Multi-chiral doublets in one single nucleus 
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Abstract 
Adiabatic and configuration-fixed constraint triaxial relativistic mean field (RMF) approaches 
are developed for the first time and a new phenomenon, the existence of multi-chiral doublets 
(MxD), i.e., more than one pairs of chiral doublets bands in one single nucleus, is suggested for 
106Rh based on the triaxial deformations together with their corresponding proton and neutron 


configurations. 
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The handedness or chirality is a subject of general interests in molecular physics, ele- 
mentary particles, and optical physics. The occurrence of chirality in nuclear physics was 
suggested in 1997 [1] and the predicted patterns of spectra exhibiting chirality, i.e., the chiral 
doublets bands, were experimentally observed in 2001 [2]. 

Since the pioneer work on chirality in nuclear physics, lots of efforts have been made to 
understand the new phenomena and explore their possible existence in the nuclear chart, 


e.g., [3, 4, 5, 6, 7, 8, 9, 10]. Experimentally, the chiral doublets bands have been identified in 


—1 


many nuclei in A ~ 130 mass region with the earlier suggested configuration th 1/2 & vh, Jo 


A ~ 100 with 795 & vhi1j2, and A ~ 190 with thg/2 & Vis. On theoretical aspect, the 
chiral symmetry breaking was firstly predicted in the particle-rotor model (PRM) and tilted 
axis cranking (TAC) approach for triaxially deformed nuclei |1]. It has been investigated later 
in hybrid Woods-Saxon and Nilsson model combined with shell correction method [7] as well 
as the Skyrme-Hartree-Fock cranking approach [8], and its selection rules for electromagnetic 
transitions have been discussed in a simple PRM [9]. 

For triaxially deformed rotational nucleus, the collective angular momentum favors align- 
ment with the intermediate axis, which in this case has the largest moment of inertia. Mean- 
while, the valence particle and hole angular momentum vectors align along the nuclear short 
and long axis, respectively. These orientations maximize the overlap of the particle densities 
with the triaxial core and minimize the interaction energy. The three mutually perpendic- 
ular angular momenta can be arranged to form two systems with opposite chirality, a left- 
and a right-handedness. They are transformed into each other by the chiral operator which 
combines time reversal and spacial rotation of 180°, x = 7R (n). The breaking of the chiral 
symmetry in atomic nucleus is observed due to the quantum tunnelling between the systems 
with opposite chirality. 

'The description for the quantum tunnelling of the chiral partners is beyond the mean field, 
as the usual cranking approach is a semiclassical model and the total angular momentum is 
not a good quantum number [11]. In contrast, the PRM is better suited for this purpose, 
which however is so far confined for one particle and one hole configurations only [1, 9, 10]. 
Its generalization for more particles and/or holes are still under development. Furthermore 
the deformation y and configurations in PRM are not self-consistent rather as inputs of the 
model. 


The relativistic mean field (RMF) theory has received wide attention due to its success in 


describing the properties of nuclei and many nuclear phenomena for the past years [12, 13]. 
It is interesting to search for nuclei with triaxial deformation and configurations with not 
only one particle and one hole but also multi-particle and multi-hole suitable for chirality in 
RMF theory. A multi-dimensional microscopic cranking RMF model is very time consuming 
and has only been applied in magnetic rotation so far [14]. 

In this Letter, we will develop the adiabatic and configuration-fixed constraint triaxial 
RMF approaches to investigate the triaxial shape coexistence and the possible chiral doublets 
bands in A ~ 100 mass region. The existence of multi-chiral doublets (MyD), i.e., more than 
one pairs of chiral doublets bands in one single nucleus, will be suggested via the examining 
the deformation and the corresponding configurations. 

In RMF theory, the nuclei are characterized by an attactive scalar field S(r) and a 


repulsive vector field V(r) in the Dirac equation, 
{-ia-V+V(r)+6[M + S(r) Hi = Eii, (1) 


which can be solved by expanding separately the upper and lower components of the spinor 
ib; in terms of eigenfunctions of the three-dimensional deformed oscillator in Cartesian co- 


ordinates ¢,(r) and its time reversal state ¢3(r) = T(r) with T = iG, K , 
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In order to avoid the complex matrix diagonaliztion problem in Ref. [15], a(r) is written 


as a product of three Hermite polynomials [16], 
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The meson fields are expanded similarly with the same deformation ĝo and ^ for the basis, 
but an oscillator length smaller by a factor of v2 than that for the nucleons (bg = bo/V2) 
has been used in order to simplify the calculations and to avoid the necessity of additional 
parameters. The oscillator frequency is hwy = 41A-!/? MeV. 

In order to check the convergence of the results with the number of expanded oscillator 


shells for fermions nor, the binding energies (upper) and the deformations 6 (middle) and 


^ (lower) calculated with effective interaction PK1 [17] for 1°°Rh as functions of no; are 
presented in Fig. 1. The legends (filled circles, open circles, squares, triangle ups, triangle 
downs, and diamonds) represent the results obtained in basis with different deformation (9— 
0.0, 0.1, 0.2, 0.3, 0.4 and 0.5, respectively). It shows that as long as nor > 10, the binding 
energies as well as the deformations 6 and y obtained are almost the same. Furthermore 
they are independent on the deformation 6o of the basis. Similar convergence check for the 
bosons has also been done. In the following, a spherical basis with 12 major oscillator shells 
for fermions and 10 shells for bosons will be used, which gives an error less than 0.196 for 
the binding energy. 

In general the triaxial RMF calculation leads to only some local minima. In order to get 
the ground state for triaxial deformed nucleus, constraint calculations are necessary and in 
principal such calculations should be carried out in two-dimensional 8 and y plane. As such 
two-dimensional constraint calculations turn out to be expensive even for modern computer 
facilities, therefore the constraint calculations with (o + 203), i.e., 8?, are carried out 
to search for the ground state for triaxially deformed nucleus. 

The energy surface and the deformation y as functions of deformation f in adiabatic con- 
straint triaxial RMF calculation with PK1 for Rh are presented as open circles in Fig. 2(a) 
and (b) respectively. There are some irregularities in the energy surface. Furthermore some 
local minima are too obscure to be recognized and it is technically difficult to understand 
their corresponding single particle configurations. We therefore performed the configuration- 
fixed constraint calculation similar to what have been done in the non-relativistic case [18]. 
Starting from any point in the energy surface in the adiabatic constraint calculations, the 
configuration-fixed constraint calculation requires that the single particle orbits occupied 


are fixed during the constraint calculation, i.e., 


(lE + 9B) oi (8))] e 1. (4) 


The energy surfaces and the deformations y in configuration-fixed calculations with PK1 for 
106Rh are given as solid lines in Fig. 2(a) and (b) respectively. For each fixed configuration, 
the constraint calculation gives a continuous and smooth curve for the energy surface and 
the deformation ^ as a function of deformation 8. The irregularities in the adiabatic energy 
surface disappear. In comparison, the minima in the energy surfaces of the configuration- 


fixed constraint calculations become obvious, which are respectively represented by stars 
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and labelled as A, B, C, D, E, F and G. Their corresponding deformations 8 and y together 
with their binding energies are respectively given in Fig. 2(a) and (b). It is interesting to 
note that, for each fixed configuration, the deformation 7 is approximately a constant (as in 
Fig 2(b)), which means that the deformation y is mainly determined by its corresponding 
configuration. 

The energies for these minima including the ground state are within 1.3 MeV to each 
other but correspond to different deformations 6 and y, which is a good example of the 
shape coexistence. The shape coexistence here is different from the spherical, oblate and 
prolate shape coexistence, e.g., in neutron-deficient Pt, Hg and Pb isotopes. It is the triaxial 
shape coexistence, i.e., for the ground state A: the binding energy E= 903.92 MeV (the data 
906.72 MeV), deformation 6 = 0.27 and y = 24.7 °, the excited minima B: E— 903.82 MeV, 
B = 0.25 and. y = 23:3 °, C: E= 903.28 MeV, 6 = 0.30 and y = 22.0 °, and D: E= 902.69 
MeV, 8 = 0.22, y = 30.8°. All the states A, B, C, and D have deformation @ and y 
suitable for chirality |1, 10]. As these states are all in one single nucleus, if chiral doublets 
bands can be built on these states, it may lead to a new phenomenon, the existence of 
MyD in one single nucleus. Therefore, in the following, we will investigate their proton and 
neutron configurations in detail to see whether the particle and hole configurations required 
by chirality are available. 

Performing the configuration-fixed constraint calculations for the ground state, the single 
particle levels can be obtained as functions of deformation 6. The difference in single 
particle levels obtained by choosing other minima is negligible. The neutron and proton 
single particle levels obtained in such a way are presented in Fig. 3. The positive (negative) 
parity states are marked by solid (dashed) lines and the occupations corresponding to the 
minima in Fig. 2 are represented by filled circles (two particles) and stars (one particle). 
The corresponding quantum numbers for spherical case are labelled at the left side of the 
levels. As the energy curves become very stiff for small deformation in Fig. 2, the present 
configuration-fixed constraint calculations cannot be performed for deformation 6 < 0.06. 
Therefore we cannot distinguish the occupation of some low-7 orbits, e.g., the last occupied 
neutron orbit for state A may come from 2d5/2 or 2d3/2, as marked in the figures. However, 
as the nuclear chirality is essentially determined by the high-7 orbits, such illegibility does 
not influence the conclusions here. 


The binding energies, deformations 6 and y as well as the corresponding configurations 


extracted from Fig. 3 for the minima A, B, C, and D in !°°Rh obtained in the configuration- 
fixed constraint triaxial RMF calculations with PK1 are listed in Table. I. Except the state 
D, in which there is no high-7 neutron valence particle, we found the high-7 proton and 
neutron configuration for the ground state (state A) as 7(19goj3) ? & v(1hi1/5)^, state B as 
m(1g9/2) ? &v(1hi1/5)!, and state C as z(1ggj3) ? &v(1h113)?, respectively. All of them have 
high-7 proton holes and high-j neutron particles configurations, which together with their 
triaxial deformations favor the construction of the chiral doublets bands. It is interesting 
to note that the states A and B compete strongly with each other in energy. However, due 
to different parities, the states A and B do not mix up and are possible to produce the new 
phenomenon MyD. So far, a pair of chiral doublets bands with negative parity have been 
observed in !??Rh [5]. It will be interesting to search for other chiral doublets bands in this 
nucleus. 

Similar to Rh, detailed constraint calculations are also performed for other isotopes 
in A ~ 100 mass region and the possibilities of MyD exist in other nuclei as well. The 
results will be published elsewhere in details. Here in Fig. 4, the deformations 6 and ^ for 
ground states in 95-114 Rp, 1027116 Ag and 1007115] isotopes are presented. The shaded area 
represents the favorable deformation ^ for nuclear chirality. The nuclei, '^Rh and '°Rh, 
in which the chiral doublets bands have been observed [4, 5], are marked as filled circles. 
Apart from “Rh and Rh, the favorable deformation y for chirality has been found in 
102,108,110Rh, 108-112 Ag and ‘In, which implicates that more chiral doublets bands can be 
expected in A ~ 100 mass region. For each isotope chain, the deformation y varies with the 
neutron number due to different occupation in v/11/5. 

In summary, adiabatic and configuration-fixed constraint triaxial RMF approaches are 
developed for the first time to investigate the triaxial shape coexistence and possible chiral 
doublets bands. A new phenomenon, the existence of multi-chiral doublets (MyD), i.e., 
more than one pairs of chiral doublets bands in one single nucleus, is suggested in Rh 
by examining the deformation and the corresponding configurations. Similar investigations 
have also been done for Rh, Ag and In isotopes and also with other effective interactions. 


The observation of MxD is very promising in this mass region. 
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TABLE I: The binding energies, deformations 8 and y and the corresponding configurations for 


the minima A, B, C, and D in !°°Rh obtained in the configuration-fixed constraint triaxial RMF 


calculations with PK1. 


903.9150 7(199/3) ? & v (111,3)? [2ds/3)! or (2d3/5)!]] 
r3 903.8196 7199/3) ? & v { (Lhy1/2)"[(2d5/2)? or (2d3/5)] j 


C 903.2790 7(19gj3) ? & v(1h113)? 
902.6960 7 (19/2) ? & v[(2ds3)? or (2da/5)] 
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FIG. 1: Total binding energy, deformation f and y calculated with PK1 for !°°Rh as functions of 
the number of expanded oscillator shells for fermions nor. The filled circles, open circles, squares, 
triangle ups, triangle downs, and diamonds represent the results as the deformation of basis 6o 


equals to 0.0, 0.1, 0.2, 0.3, 0.4 and 0.5, respectively. 
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FIG. 2: (color online) The energy surfaces (a) and the deformations y (b) as functions of defor- 
mation f in adiabatic (open circles) and configuration-fixed (solid lines) constraint triaxial RMF 
calculation with PK1 for !°°Rh. The minima in the energy surfaces are represented as stars and 
labelled respectively as A, B, C, D, E, F and G. Their corresponding deformations @ and y together 


with their energies are respectively given in (a) and (b). 
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FIG. 3: (color online) The neutron and proton single particle levels obtained in configuration-fixed 
constraint triaxial RMF calculations with PK1 for !°°Rh as functions of the deformation 8. Positive 
(negative) parity states are marked by solid (dashed) lines. The occupations corresponding to the 


minima in Fig. 2 are represented by filled circles (two particles) and stars (one particle). 
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FIG. 4: (color online) The deformations £8 and y for the ground states in Rh (upper), Ag (middle) 
and In (lower) isotopes in constraint triaxial RMF calculations with PK1. The shaded area repre- 
sents the favorable deformation ^ for nuclear chirality. The nuclei !°4Rh and !°°Rh, in which the 


chiral doublets bands have been observed, are marked as filled circles. 
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